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Abstract. - We calculate the noise spectrum of the electrical current in a quantum point 
contact which is used for continuous measurements of a two-level system (qubit). We generalize 
the previous results obtained for the regime of high transport voltages (when V is much larger 
than the qubit's energy level splitting B (we put e = ?!. = 1)) to the case of arbitrary voltages 
and temperatures. When V B the background output spectrum is essentially asymmetric 
in frequency, i.e., it is no longer classical. Yet, the spectrum of the amplified signal, i.e., 
the two coherent peaks at uj = ±B is still symmetric. In the emission (negative frequency) 
part of the spectrum the coherent peak can be 8 times higher than the background pedestal. 
Alternatively, this ratio can be seen in the directly measureable excess noise. For V < B and 
T = the coherent peaks do not appear at all. We relate these results to the properties of 
linear amplifiers. 



Introduction. - For quantum information technology it is necessary to investigate prop- 
erties of real physical systems used as quantum detectors. Certain quantum algorithms require 
an efficient (single-shot) read out the final state of a qubit. This can be done by either strongly 
coupled threshold detectors (see e.g., Refs. [1,2]), or by "measurement in stages" strategy [3]. 
For weakly coupled detectors the only way to perform single-shot measurements is to be in 
the qunatum- non-demolition (QND) regime, i.e., by measuring an observable which commutes 
with the Hamiltonian and is, thus, conserved. In the solid state domain this regime has been 
investigated in, e.g., Refs. [4-6]. 

In this letter we concentrate on continuous weak non-QND measurements (monitoring) of 
the coherent oscillations of a qubit (two-level system, spin- 1/2). This regime was the main 
focus of Refs. [7-10]. It is realized, e.g., for the transverse coupling between the spin and the 
meter, e.g., when the effective magnetic field acting on the spin is along the x-axis while cr^ is 
being measured. In this case one observes the stationary state properties of the system, after 
the information about the initial state of the qubit is lost. Thus, this regime is not useful for 
quantum computation. Yet, studying the properties of the meter in the stationary monitoring 
regime, one can obtain information nessecary in order to, later, employ the meter in the QND 
regime. Another motivation for our study comes from the recent activity in the STM single 
spin detection (see, e.g., Ref. [11,12]). 
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Without monitoring and without couphng to other sources of dissipation the observable 
(Jz would show coherent (Larmor) oscillations. When subject to monitoring these oscillations 
give rise to a peak in the output spectrum of the meter at the Larmor frequency. The laws 
of quantum mechanics limit the possible height of the peak. In the case of a 100% efficient 
(quantum limited) detector and when all the noises are white on the frequency scale B the 
peak can be only 4 times higher than the background noise pedestal [7] . Inefficiency of the 
detector reduces the height of the peak further. 

Usually the analysis of the continuous measurements in voltage driven meters is limited 
to the case V ^ B [4,7,8]. The output noise spectrum in this regime is almost symmetric 
(classical) at frequencies of order and smaller than B. In this letter we remove the restriction 
V ^ B. At low voltages, V ^ B, the output noise is essentially asymmetric, i.e., the output 
signal is quantum. In other words, we have to differentiate between the absorption (w > 0) and 
the emission (w < 0) spectra of the detector (see, e.g. Ref. [13]). We show, however, that the 
qubit's contribution to the full output noise as well as to the experimentally accessible excess 
noise is symmetric. We calculate this contribution for arbitrary voltage and temperature. In 
the excess noise, which is obtained by subtracting the equilibrium detector noise {V = 0) from 
the output at y ^ 0, the peak to background ratio can reach 8 for ^ ~ B. 

General considerations. - We start from the general theory of linear amplifiers [5,6,9,14] 
which applies in the regime of weak continuous monitoring. The Hamiltonian of the whole 
system including the amplifier (meter) reads H = Hector + Hqs+caQ, where a is the measured 
observable of the small quantum system governed by H^s, Q is the input variable of the 
amplifier governed by iJmeter, and c is the coupling constant. The meter is necessarily driven 
out of equilibrium. We study the output variable of the meter /. The stationary average 
value (/) is only slightly changed by the presence of the qubit. Much more interesting is 
the spectrum of fluctuations {SI^) = J dt {5I{t)SI{0)) e^'^* . While it is convenient to discuss 
physics in terms of the symmetrized Si{ui) = (l/2)[((5/^) + {61^^)] and anti-symmetrized 
Ai{aj) = (l/2)[((5/2) — {Sl'^^)] correlators, the calculations are more convenient in terms of 
the Keldysh-time-ordered Green's functions (see, e.g., Ref. [15]). Thus we define Gi{t,t') = 
—i{TKSI{t)6I{t')). This is a 2 x 2 matrix as both t and t' can belong either to the forward 
or to the backward Keldysh contours [15]. We have the two basic components iGf{t — 
t') = {6I{t)SI{t')} and iGf{t - t') = {5I{t')6I{t)) from which all others are built. The 
retarded and advanced components are defined as G^{t — t') = 9{t' — t)[Gf{t — t') — Gf{t — 
t')] and G^{t - t') = ~6{t - t')[G^{t - t') - Gf{t - t')]. These two components describe, 
usually, the response of 7 to a perturbation coupled to I. The Keldysh component defined 
as G/ (t — t') = Gf{t — t') + Gf(t — t') is related to the symmetrized correlator. It is easy 
to obtain the following relations: iG^{uj) - iGf{ijj) = 2Ai{uj) and iGj (w) = 2Si{uj). We 



will also need GiQ{t,t') = -i{TK6I{t)SQ{t')), Gqiit^t') = -i{TKSQ{t)SI{t')), and n{t,t') = 



—i{TKSa{T)Sa{T')) . Various components of these Green's functions are defined analogously 
to those of Gj. 

We assume that one is allowed to use Wick's theorem for the operators /, Q, and a. 
Frequently, even if Wick's theorem does not apply, one can still use it for the lowest (in the 
coupling constant) calculations and show that the corrections arc of the higher order. We 
return to this subject later. Since the dynamics of the measured system changes substantially 
as a result of measurement while the meter's one is perturbed weakly, we use the full ( "thick" ) 
Green's function 11, while for Gjq and Gqi one keeps the unperturbed values. Then the lowest 
order irreducible correction to the Green's function Gi{t,t') reads 




(1) 
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In the stationary regime this gives 















GW) 



(2) 



For the Green's function Gqi we have Gqj{lj) = Gj-Q{—ijj), Gqj{uj) — GfQ{—uj), and 
G^q{—uj). We introduce the notations X{uj) — cG^Q(a;), A'(ijj) = cGqj{uj), where 



GW) 



A is the direct gain (amphfication coefficient) of the amphfier, while A' is the inverse gain. As 
X{t) and X'{t) are real, X{—lu) — X*{lu) and X'{—uj) = X'*{uj). Thus we obtain 



5Gf{uj) = X{uj)X'{uj)U^iuj) , 
SGfiuj) = -2iSSi{uj) = |A(w)p U^{uj) + 2iclm [X{lu) U^{lu) G^j{uj)] 



(3) 
(4) 



and SGfiuj) ^ {6G^\uj)\* . We have also used G§j(w) = G%{-uj) = --[Gf^iLu)]* . The first 
term of (@J corresponds to the noise of the small system "amplified" by the meter. The second 
term is needed to fulfill the fluctuations-dissipation relation at equilibrium. We will see that 
it is also important at low voltages, i.e., when the detector is not driven far enough from 
equilibrium. For good amplifiers the inverse gain vanishes. A' = 0, and we obtain 6G^ ~ 0. 
Thus the contribution to the output correlator S (SI^) = iSGf = {i/2){SGf + SGf ~ SGf) = 
{i/2)5G^ — SSi{uj) is symmetric in frequency, i.e., SAi{uj) — 0. Vanishing of A' also means 
that further amplifiers using / as an input will not add to the back-action. 

Spin's dynamics. - When an observable of a qubit (a component of spin-1/2) is being 
measured we can assume without loss of generality a = az. The spin is subject to an (effective) 
magnetic field B, i.e., Hqs = — {1/2)B a. Its dynamics is, thus, obtained from the Hamiltonian 



(5) 



where we have put c — —1/2 so that Q can be interpreted as fluctuating magnetic field. 
We exclude the case B \\ z, in which the measured observable a — az commutes with the 
Hamiltonian and, thus, is conserved. This regime is known as the quantum-non-demolition 
(QND) one and has been treated, e.g., in Refs. [4-6]. In all other cases the stationary state is 
reached after some transient period and we can study the output spectrum of the meter. 

For simplicity we assume no extra dissipation sources acting on the qubit except for the me- 
ter. The spin's Green functions (correlators) are obtained within the standard Bloch-Redfield 
approach [16, 17] which is applicable as long as the dissipation is weak (see below). Within 
this approach one, first, calculates the markovian evolution operator for the spin's density 
matrix, and, then, uses the "quantum regression theorem" [13] to obtain the correlators. For 
this lowest order perturbative (in the spin-meter coupling) calculation one only needs to know 
the (unperturbed by the spin) fluctuations spectrum, (Q^). One, then, obtains 



U^{uj) = -isin^ 



2r2 



2To 



u-Bf+Tl {u: + BY+Tl 
1 1 



I cos 



4ri 



^2 + r2 



[1 - {-sn . 



LO - B + iT2 + B + iT2 



(6) 



where 9 is the angle between B and z. The stationary spin polarization along B is given 
by (a^) = h{B), where h{io) = Aq{uj) / Sq{uj), while Sq{u) = {l/2)[{Ql) + {Q^_^)] and 
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by: 



(1/2)[(Q2> - {Q-J]- The relaxation (Fi) and the dephasing (r2) rates are given 



ri^{l/2)sinHSQ{Lu = B) , r2 = (l/2)ri + (l/2)cos2(?5Q(cj = 0) 



(7) 



The apphcabihty condition of the Bloch-Redfield approach is ri,r2, SB <C B, where SB is 
the renormahzation of the energy sphtting (Lamb shift). If one treats the Lorentzians in © 
as true delta functions, one can derive from the second equation of |(BJ a relation resembling 
the fluctuation-dissipation theorem: — n^(aj) w h{uj) {lj) . Note, that, although we 

use the "diagrammatic" language of Keldysh Green functions, Eqs. ^ are obtained without 
any diagrams or assumptions about the applicability of Wick's theorem. 

Spin's contribution to the output spectrum. - Substituting Eqs. © into Eq. Q we obtain 
spin's contribution to the symmetrized output spectrum of the meter: 



dSiiuj > 0) 



sm 



+ sin^ 



uj- B 



cos^ e ■ 



{i.-B)^ + rl 
2ri 



|AH|2-^Re[AMG^,H] 
^Im[AHG^,H] 



[l~h'{B)] |AM|^ 



(8) 



We assume that A(a;) and Gqj(ijj) are smooth near w — zLB and u; = 0. Then, the first term 
of JS]) gives two peaks near uj — zLB with width The third term gives a peak at w = 
with widt Fi. The second term of (jSJ gives the Fano shaped contributions near lu = ±B. The 
simplest situation arises when = 7t/2 and h{B) (at very high transport voltages the 
effective temperature of the spin is infinite and (cr^) — > 0). Then only the peaks at w = ±B 
survive with the height SSiiB) = |A(B)|Vr2 4\X{B)\^ / Sq{B). The peak to pedestal ratio 
SSi{B)/SiiB) = A\X{B)\^/{SQiB)Si{B)) was shown [7-10] to be limited by 4. Below we 
investigate the coherent peaks at arbitrary voltage and temperature for a specific example of 
a meter with the purpose to explore the effect of the rest of the terms in . 

Quantum Point Contact (QPC) as a meter. - The QPC devices are known to serve as 
effective meters of charge (see, e.g., Refs. [18-21]). The conductance of the QPC is controlled 
by the quantum state of a qubit. We focus on the simplest limit of a tunnel junction when 
the transmissions of all the transport channels are much smaller than unity. This model has 
previously been used by many authors [22-24]. It also corresponds to the model considered 
in Ref. [12] in the regime of lead electrons fully polarized along the z axis, niR = niL = z. 
The tunnel junction limit is described by the following Hamiltonian 



E(^o + ha,){cl c; + h.c.) - (1/2) 



(9) 



The transmission amplitudes t^ and ti are assumed to be real positive and small (tunnel 
junction limit). We also assume ii ^ to to be in the linear amplifier regime. 

For brevity we introduce the operator X = ^ c\ c^ and then j = i{X — X^) and q = 
{X ~\-X'^). The current operator is given by / = [t^ +tiaz)j, while the tunneling Hamiltonian 
is Ht — (to + tiaz)q. We see that the analysis of the amplifiers presented above cannot be 
directly applied. First, the interaction term between the spin and the QPC, i.e., ti qaz (thus, 
in our case Q = —2tiq), is not the full interaction vertex but rather a part of H^. Second, the 
current operator / contains the spin's operator cjj explicitly. One possible way to resolve these 
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a) 



c) 




Fig. 1 - The diagrams leading to Eq. I|13^ . The loops represent the Green functions Gjq and Gqj. 
Each of these Green's functions is actually a combination of two electronic ones (solid lines). The 
dashed line is the spin's Green function. 



difficulties (see e.g., [7,9]) is to include the spin-independent part of H^, namely to q, into the 
zeroth-order Hamiltonian. This amounts to working in the basis of scattering states. Here we 
adopt a simpler procedure suitable for QPC's in the tunneling regime. We expand in the full 
Ht and keep all the terms up to the order t\ . For this we need the following zeroth-order 
Green's functions: G,, = -i{TKq{t)q{t')) , G,, = -i{TKjit)j{t')) , and G,, = -i{TKjit)q{t')), 
Gqj = —i{TKq{t)j{t')). For w <C -D, where D is the electronic bandwidth (the Fermi energy) 
we obtain 

LU + i... 2s(lu) 



GM = G,,H = - j , (10) 

G,q{u;)-~Gq,[u;)-rj^ -V{1 - tO{u;/D)) ) ' 

where 77 = 27rpLPR- We have also introduced the two following functions: 

V + uj ^ V + Lu V-oj ^ V-uj 
s/a (w) EE coth ± coth . (12) 

In Eq. Hl()(l ... stand for the real part of the retarded (advanced) components. The factors 
1 ± iO{uj/D) in are responsible for making the functions G^^{t) and G^^{t) causal. As 
we are interested in the low frequencies [to ^ D) we approximate those factors by 1. 

Peaks in the output noise spectrum. - We combine the Green functions Gjq{uj), Gqj{uj), 
and n(cL;) (see Eq.^J into diagrams to calculate the qubit's contribution to the current-current 
Green function SGi. The Wick theorem does not apply to the spin operators. However, using 
the Majorana representation of the spin operators and recently proved useful identities (see 
Refs. [25,26]) we are able to show that in order t^tl the answer is given by the diagrams shown 
in Fig. HI We obtain SGi = t^tl [Gjq{uj) + Gfg{0) ■ i] H(w) [Gqj{uj) + Gfj{0) ■ i] , which can 
be rewritten as 

xr^ M I 2^ f V a{uj) \ / H« HK \ / ~a{oj) \ 

5Gi^{il^')9m[ V )[ i) HA j V ) ^ (1^) 

where we have defined the conductances as = 2'Kritf^ and gi = 2T:rit\ It is worth comparing 
Eqs. and (jSJ. Even though the simple formalism leading to Eq. ^ was not directly 
applicable in our case, the result looks very similar. We can interpret, therefore, A = 
(I/tt) y^ogTF, A' — 0, and cG^q(w) — (I/tt) ^goffi a{u}). Thus the tunnel barrier possesses 

R / A 

the property A' = at all frequencies. This assures that SG^ = and the contribution to 
the current-current correlator is symmetric: 



5S'/(w>0) = (I/tt^) 30.91 ^^ sin^^* 



{uj - Bf + r 



2 



a{uj)h{B) 
V 



2r 

+ (l/^2)^^^^^2^^g2^_l [l_/,2(5)] . (14) 
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Fig. 2 - Output noise power for 6 = n/2, go = 0.25, gi = 0.05, T = 0.05B, and V = 1.3B. Solid line: 
the full correlator; dashed line: the excess spectrum. 

Fig. 3 - Same parameters but 9 — 2tt/3. Only the full correlator shown. 

We note that in our example Sq{lj) = 2itiG^g{uj) = {2/n) gis{uj) and Aq{lj) = (2/11) giLO. 
Then we obtain h{B) = B/s{B), Ti = (I/tt) gi sin^ d s{B), and T2 = (l/27r) gi sin^ 6 s{B) + 
(I/tt) gi cos^ 6 s(0). Note, that no Fano shaped contributions appear due to the fact that both 
A(a;) and G^q{lu) are real. The Lorentzians in Eq. 1)14(1 coincide with the ones obtained in 
Refs. [7-10]. The new result is the reduction factor for the peaks at w = zkB in the square 
brackets. This factor simplifies for T = 0. Then, iiV>B, it is given by (1 — B^/V^), while 
for F < i? it is equal to 0. In the last case the measuring device can not provide enough 
energy to excite the qubit and, therefore, the qubit remains in the ground state and does not 
produce any additional noise (see also Ref. [12]). The ratio between the peak's height and 
the pedestal's height is different for positive and negative frequencies. In the limit gi <C go 
we obtain (SI^) « itlG>Juj) = (l/27r) gQ{s{uj) + uj). Thus, for T = 0, and B <V we obtain 



(SIIb) ~ (l/27r) 30^(1 ± B/V) and S(Sllg) = SSi{±B) « (2/7r) .goT^(l - B^ /V^) and 



For B ^ V the ratio for the negative frequency peak reaches 8. In this limit, however, the 
peak's hight is zero. For symmetrized spectra the maximal possible ratio is 4 (Ref. [8]). An 
interesting question is what exactly can be observed experimentally. If the further detection 
of the output noise is passive, like the photon counting in fluorescence experiments, one can 
only measure what the system emits, i.e. the noise at negative frequencies [27-29]. Moreover, 
in our example, the excess noise, i.e., {6Ij.g){V) — {6Ij.g){V = 0), is symmetric. As shown 
in Ref. [30], if the excess noise is symmetric, it can be effectively measured even by a finite 
temperature LCR filter. In Figs. 1^151 we show examples of output noise spectrum and of the 
corresponding excess noise spectrum. 

Conclusions. - We have calculated the output noise of the point contact used as a 
quantum detector of qubit's coherent oscillations for arbitrary voltage and temperature. In 
the regime eV ~ B and T ^ B the output noise is essentially asymmetric. Yet, the qubit's 
oscillations produce two symmetric peaks at a; = ±B and also a peak at w = 0. Due to the 
vanishing of the inverse gain (A' = 0) the peaks at oj — ±B have equal height and, therefore, 
the negative frequency peak is much higher relative to it's pedestal than the positive frequency 
one. The peak/pedestal ratio can reach 8. This can be observed by further passive detectors, 
which measure what the system emits, or by measuring the excess noise. The results of this 
paper are obtained for the simplest and somewhat artificial model of a quantum detector, a 




(15) 
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QPC in the tunneling regime. It would be interesting to perform analogous calculations for 
more realistic detectors like SET's or QPC's with open channels (see e.g., Refs. [5,29,31]). 

* * * 

We thank Yu. Makhlin, Y. Levinson, and L. Bulaevskii for numerous fruitful discussions. 
A.S. was supported by the EU 1ST Project SQUBIT and by the CFN (DFG). D.M. and I.M. 
were supported by the U.S. DOE. 

REFERENCES 

[1] D. Vion, A. Aassimo, A. Cottet, P. Joyez, H. Pothier, C. Urbina, D. Esteve, and M. H. Devoret, 
Science 296, 886 (2002). 

[2] I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, Science 299, 1869 (2003). 
[3] O. Astafiev, Yu. A. Pashkin, T. Yamamoto, Y. Nakamura, and J. S. Tsai, cond-mat/ 04026 19 

(2004). 

[4] A. Shnirman and G. Schon, Phys. Rev. B 57, 15400 (1998). 

[5] A. A. Clerk, S.M. Girvin, and A.D. Stone, Phys. Rev. B 67, 165324 (2003). 

[6] M. H. Devoret and R. J. Schoelkopf, Nature 406, 1039 (2000). 

[7] A. N. Korotkov and D. V. Averin, Phys. Rev. B 64, 165310 (2001). 

[8] A. N. Korotkov, Phys. Rev. B 63, 085312 (2001). 

[9] D. V. Averin, In J. R. Friedman and S. Han, editors. Exploring the Quantum-Classical Frontier 

Commack, NY; cond-mat/0004364 (2002). Nova Science Publishers. 
[10] R. Ruskov and A. Korotkov, Phys. Rev. B 67, 075303 (2003). 
[11] A. V. Balatsky and 1. Martin, cond-mat/0112407 (2001). 
[12] L. N. Bulaevskii, M. Hruska, and G. Ortiz, Phys. Rev. B 68, 125415 (2003). 
[13] C. W. Gardiner and P. Zollcr, Quantum Noise, Springer 2 edition (2000). 
[14] V. B. Braginsky and F. Ya. Khalili, Quantum measurement, Cambridge University Press 

Cambridge (1992). 
[15] J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
[16] F. Bloch, Phys. Rev. 105, 1206 (1957). 
[17] A. G. Redfield, IBM J. Res. Dev. 1, 19 (1957). 

[18] M. Field, C. G. Smith, M. Pepper, D. A. Ritchie, J. E. F. Frost, G. A. C. Jones, and D. G. 

Hasko, Phys. Rev. Lett. 70, 1311 (1993). 
[19] D. Sprinzak, Yang Ji, M. Heiblum, D. Mahalu, and H. Shtrikman, Phys. Rev. Lett. 88, 176805 

(2002). 

[20] E. Buks, R. Schuster, M. Heiblum, D. Mahalu, and V. Umansky, Nature 391, 871 (1998). 
[21] J. M. Elzcrman, R. Hanson, ,J. S. Groidanus, L. H. Willcms van Bovcrcn, S. Dc Franceschi, 

S. Tarucha L. M. K. Vaudcrsypou, aud L. P. Kouwenhoven, cond-mat/0212489 (2002). 
[22] S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 
[23] A. N. Korotkov, Phys. Rev. B 60, 5737 (1999). 
[24] H. S. Goan and G. J. Milburn, Phys. Rev. B 64, 235307 (2001). 

[25] W. Mao, P. Coleman, C. Hoolcy and D. Langrcth, Phys. Rev. Lett. 91, 207203 (2003). 

[26] A. Shnirman and Yu. Makhlin, Phys. Rev. Lett. 91, 207204 (2003). 

[27] G. B. Lesovik and R. Loosen, JETP Lett. 65, 295 (1997). 

[28] U. Gavish, Y. Levinson, and Y. Imry Phys. Rev. B 62, R10637 (2000). 

[29] R. Aguado and L. P. Kouwenhoven, Phys. Rev. Lett. 84, 1986 (2000). 

[30] U. Gavish, Y. Imry, Y. Levinson, and B. Yurkc, In Y.V. Nazarov and Y. Blanter, editors. 

Quantum Noise in Mesoscopic Physics. Kluwer (2003). 
[31] S. Pilgram and M. Biittiker, Phys. Rev. Lett. 89, 200401 (2002). 



